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Abstract 

We propose a scheme for realizing the Ising spin-spin interaction and atomic cluster states 
utilizing trapped atoms in coupled microcavities. It is shown that the atoms can interact with each 
other via the exchange of virtual photons of the cavities. Through suitably tuning the parameters, 
an effective Ising spin-spin interaction can be generated in this optical system, which is used to 
produce the cluster states. This scheme does not need the preparation of initial states of atoms 
and cavity modes, and is insensitive to cavity decay. 
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Strongly correlated many-particle systems have been extensively explored in condensed 



matter physics, cold atoms 
lattices [sl, 



l|, |2|, and recently in optical system such as coupled microcavity 
Compared to other strongly correlated many particle systems, 



an optical system has the advantage of easily addressing individual lattice sites with optical 
lasers. Because of the size and separation of the microcavities, arbitrary lattice geometries 
can be arranged in this system. Therefore, it offers the ability to experimentally observe 
quantum-many-body phenomena and realize quantum information processing. Various tech- 
nologies have been employed in this field, including fiber coupled toroidal microcavities [lO|, 
arrays of defects in photonic band gap materials and superconducting qubits coupled 
through microwave stripline resonators (l2[]. 

On the other hand, interacting qubits are particularly important in quantum information 
processing. Lattices of interacting spins or qubits can be utilized to generate highly entangled 



states, such as cluster states 



13 



14| . It has been shown that this class of entangled states are 



much more entangled than the GHZ states, thus have a high persistence of entanglement 13|, 



14l |. The experimental demonstration of the violation of Bell's inequality for cluster states 



has been reported 151] • Moreover, cluster states together with local measurements are the 



resource for one-way quantum computation 



16 



171]. In the context of cavity QEDjlS 



201 ]. several proposals have been presented for producing the atomic cluster states 




However, the scalable implementation in many-particle case is difficulty. Recently in Ref. 
a very novel scheme has been proposed to simulate the dynamics of an effective anisotropic 
Heisenberg spin-1/2 chains and generate atomic cluster states in coupled cavities. However, 
there is no report of the realization of Ising interaction and cluster states in a simple enough 
atom-cavity system that does not require the complicated atomic level structure and driving 
laser configurations. 

In this paper we present a scheme for realizing the Ising spin-spin interaction and atomic 
cluster states utilizing trapped two-level atoms in coupled microcavities arranged in an array. 
We show that the atoms can interact with each other via the exchange of virtual photons of 
the cavities. Through suitably tuning the parameters, an effective Ising spin-spin interaction 
can be generated in this optical system. We discuss how to use this Ising interaction to 
produce the cluster states of atoms. This scheme does not need the preparation of initial 
states of atoms and cavity modes, and is insensitive to cavity decay. With presently available 
experimental setups in cavity QED, the implementation of this scheme is feasible. 
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As sketched in Fig. 1, an array of microcavities are coupled via the exchange of photons 
with identical two-level atom in each cavity. The ground state of each atom is labeled as 
\gj)-, and the excited state as |ej), where the index j counts the cavities. The cavity mode 
(frequency uJc) couples to the transition \g) ^ |e) (transition frequency c^o) with coupling 
constants g. Furthermore, an external strong classical field (frequency ujl) drives the same 



transition with Rabi frequencies f2[23|]. For simplicity, we only consider the one-dimensional 
array. One can generalize to the higher dimensions straightforwardly. In the rotating-wave 
approximation, the associated Hamiltonian reads (let h= 1) 

H = Ha + Hc + Hac, 

N 

Ha = ^UJoSfSj- 

i=i 

N N 

He = ujc'^ OjCij + Jc ^(aja^+i + fljaj+i), 
i=i i=i 

TV 

Hac = J2[g{a]S-+a,Sl) + n{S+e-'-^' + S-e'-^% 
i=i 

(1) 

where = \ej){gj\, Sj' = {Sj^Y, aj is annihilation operator for the photon in cavity j, and 
Jc is the hopping rate of photons between neighboring cavities. In the complete Hamiltonian, 
Ha describes the free Hamiltonian for atoms; He describes the Hamiltonian for photons in 
each cavity modes, with photon hopping between neighboring cavities; finally Hac is the 
Hamiltonian that describes the interaction between atoms and the cavities as well as the 
strong driving by the classical fields. We consider the periodic boundary conditions, then 
He can be diagonalized through the Fourier transformation. For convenience we introduce 
the notation J = [uj, 0, 0) to denote the position of the jth site where u is the length of 
the one-dimensional crystal cell and in the following chosen as unit for simplicity. Then we 
obtain He = J2k ^kO^k'^k, where tUfc = tUc + 2 Jc cos k. 

The Hamiltonian of Eq. ([T]) can be changed to a reference frame rotating with the driving 
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field frequency ujl, 

N N 

j=l k j=l 

N 

+ E + 9 St E «^^"^''']' (2) 

j=l k k 

wliere A = ujq — ujl, and 6k = ujl — Uk- In tlie following we assume uq = ujl for simplicity. 
We now switch to a new atomic basis | Ij) = + l^j)) and | tj) = '^(Idj) " then 

can rewrite H as 

AT 



k j=i 
j=l k 



•2-^ 2^ 2 

fc 



where (t| = | \ - \ Ti)(Ti I; cr+ = | Tj)(ij I; and = (cTj'y. In the interaction 

picture with respect to Hq = —J2k^kol.ak + XljLi^^l' "^^ have the following interaction 



Hamiltonian 23 1 



^11 1 



'2 2 ^ 2 



2^2^ 2 

(4) 

In the strong driving regime Q ^ (7,5^ (for all k), we can realize a rotating- wave approxi- 
mation and neglect the fast oscillating terms. Then Hj reduces to 

^ 1 

j=l k k 

(5) 

To further reduce the model, we assume 6k ^ g (for all k). Then there is no energy 
exchange between the atomic system and the cavities. We can adiabatically eliminate the 



photons from the description 



24i \. We consider the terms up to second order in the effective 
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Hamiltonian and drop the fast oscillating terms. Then we obtain the following effective 
Hamiltonian describing the Ising type spin-spin interaction in the optical system 

N 



„2„ife 



where Jz = The parameters Jz can be tuned by varying coupling strength g and 

detuning 6k- The evolution operator for the system is given by 

U{t) = e-'"'' = e-'[^f=i ^--M+il*. (7) 

It has been shown that Ising interaction can be utilized to generate a new class of multipartite 
entanglement, the so called cluster states[l^. In the following discussions we will use this 
Hamiltonian to produce the cluster states. 

In order to reduce notations we will not use the Ising interaction in Eq. but rather 
the phase gate 

^ 1 _|_ 2 1 _ z 

U,it) = e-^^'^^ with g; = -^J.t (8) 



as the elementary two-qubit interaction between neighboring atoms [13|, [IJ]. The equivalence 
between the quantum Ising interaction Hi in Eq. ([6]) and H'^ in Eq. ([8]) can be seen from 
the following discussion. From 



N 



H'j= -J^Jzt 



i=i 

N 



- W + - 0^1+1 - o^l^l+i) (9) 
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we find 



Therefore, the phase gate corresponds to the Ising interaction up to some additional 

n 

rotations around the z-axes at each qubit[14l]. The entanglement properties are determined 
by the pure Ising interaction, which are not changed by the z-rotations. 

We first consider the case of two coupled cavities. We denote = ^(| [j) + | = 
\gj), and |— j) = ^(| ij) — | Tj)) = = 1)2) for atomic states. Assume that initially the 



atom in each cavity is prepared in state = \gj). After an interaction time of t = vr/J^, 
we find that 

U'/{n/J,) = Pl_ ®f + ® al, (11) 
where P^_t = — 1) 2). Then the state evolution of the system is given by 

U,i^/J.)\+l)\+2) = ^(1 il)|-2) + I Tl)|+2)) 

= ^(Ul)^2 + lTl))(U2) + |T2)). (12) 

This state is a maximally entangled state. Up to a local unitary transformation on qubit 2, 
we can write it in the standard form. 

We now turn to the case of many coupled microcavities. Initially all the atoms are 
prepared in states = \gj){j = 1,2, ...,N). We choose t = n/ J^, then the produced state 
can be written in the following compact formjis], [l^ 

1 ^ 

l^^) = W®(U:'>l+i + lt.))- (13) 



2^ 



Then for = 3, 4, one can obtain 



1^3) = ^(|il)|i2)|i3) + |Tl)|T2)|T3))Ln. 

1^4) = ^(|il)|i2)|i3)U4) + Ul)U2)|T3)|T4) 

+ |Tl)|T2)|i3)|i4)-|Tl)|T2)|T3)|T4))L.. 



(14) 



Where "l.u." indicates the equality holds up to a local unitary transformation on one or 

nn 

more of the qubits|13l. Il4j|. {ip^) corresponds to a GHZ state of three qubits, but lip^) is not 
equivalent to a 4-qubit GHZ state. The state {ipi) is the 4-qubit cluster state. 

It is necessary to verify the approximations by numerics. We numerically simulate the 
dynamics generated by the full Hamiltonian H and compare it with the results generated by 
the effective model ([6]). As an example here we only consider the case of two atoms in two 
cavities. Initially the atoms stay in the state \g)i\g)2, and the cavity mode in the vacuum. 
In Fig. 2(a) we display the occupation probabihty Pgig2 of system in the state \gi)\g2), 
and the occupation of the photon number occupation p^Vi = {(iiO-i) calculated both from 
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the full Hamiltonian and the effective model. Fig. 2(b) shows the von Neumann entropy 
of the reduced density matrix of one spin, i.e., Ey^ = — Tr(pi logg pi(t)). A maximally 
entangled state for the atoms occurs for t = {2n + 1)tt/4Jz for any integer n. We choose 
the parameters as Q = 50 GHZ, ^ = 0.1 GHZ, = 0.02 GHZ and u^-lul = 1 GHZ. It 
can be seen that the effective model can describe the dynamics very well provided that the 
parameters are appropriately chosen. The occupations of the photon number are always 
smaller than 0.01. Discrepancies between the numerical results for the full Hamiltonian 
and the effective Hamiltonian are due to the higher order terms for the detunings and Rabi 
frequencies. However, these discrepancies are below 2% with respect to the results from the 
full Hamiltonian. 

We now consider some experimental matters. For experimental implementation, the 
parameter of the effective Hamiltonian must be much larger than decay rates of cavity 
and atomic excited states. The cavity decay has neglectable effect on this scheme. We have 
only to consider the effect of the decay of atomic excited states. This proposal requires 
that the life time of the atomic excited states be longer than the time needed to complete 
the whole procedure. For potential atomic system, Rydberg atoms are good candidate. 
Promising candidates for microcavities are photonic bandgap cavities ll|] , and toroidal or 
spherical microcavities coupled via tapered optical fibers|10|. 

In summary, we have proposed a scheme for realizing the Ising spin-spin interaction and 
atomic cluster states utilizing trapped atoms in coupled microcavities. It is shown that 
the atoms can interact with each other via the exchange of virtual photons of the cavities. 
Through suitably tuning the parameters, an effective Ising spin-spin interaction can be 
generated in this optical system. We discuss how to use this Ising interaction to produce the 
cluster states of atoms. This scheme does not need the preparation of initial states of atoms 
and cavity modes, and is insensitive to cavity decay. With presently available experimental 
setups in cavity QED, it may be implemented. 
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FIG. 1: An array of microcavities as described in the scheme. Photon hopping occurs because of 
the overlap of cavity modes of adjacent resonators. Two-level atom in each cavity is driven by 
external strong field. 
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FIG. 2: (a) The occupation probability Pgiga of the system in the state |5i)|52)(soUd hne represents 
results from the effective model, and dash line represents results from the full Hamiltonian) , and 
the photon number occupation jojv^ = (a|ai). (b) The von Neumann entropy E^iq of the reduced 
density matrix of 1 effective spin. Parameters are chosen as = 50 GHZ, g = 0.1 GHZ, Jc = 0.02 
GHZ and lOc-ujl = 1 GHZ. 
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